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Pressure can tune material's electronic properties and control its quantum state, making some systems 
present disconnected superconducting region as observed in iron chalcogenides and heavy fermion 
CeCu 2 Si 2 . For CaC 6 superconductor ( T c of 1 1 .5 K), applying pressure first T c increases and then suppresses 
and the superconductivity of this compound is eventually disappeared at about 18 GPa. Here, we report a 
theoretical finding of the re-emergence of superconductivity in heavily compressed CaC 6 . The predicted 
phase III (space group Pmmn) with formation of carbon nanofoam is found to be stable at wide pressure 
range with a T c up to 14.7 K at 78 GPa. Diamond-like carbon structure is adhered to the phase IV (Cmcm) 
for compressed CaC 6 after 126 GPa, which has bad metallic behavior, indicating again departure from 
superconductivity. Re-emerged superconductivity in compressed CaC 6 paves a new way to design new- type 
superconductor by inserting metal into nanoporous host lattice. 

The metal doped organic compounds, graphite, and solids with Tt-electrons networks, covering intercalation of 
alkali metals into buckminsterfullerene solids, metal intercalated graphite intercalation compounds (GICs), and 
alkali-metal-doped hydrocarbon materials, have recently received much attention in unveiling the new physical 
properties 1 , in particular superconductivity. So far, the highest T c achieved in this class of materials is 38 K in 
pressurized Cs 3 C 60 2-3 . In the case of GICs, superconductivity was first observed in KC S with a transition temperature 
of 0.15 K 4 . When using calcium atoms to replace the potassium atoms, which brought up to T c of 11.5 K 5,6 , making 
CaC 6 considered to be a template for exploring physical properties and superconducting mechanism of GlCs-based 
superconductors as well as improving T c value. Pressure, as a powerful tool, is often used to elevate superconductors' 
T c . Gauzzi et al concluded that a large linear increase of T c , from the ambient pressure value 11.5 K up to 15.1 K at 
7.5 GPa, followed by a sudden drop to 5 K at about 8 GPa due to a first-order structural phase transition 7 . 
Theoretically, pressured-induced structural phase transition is ascribed to softening phonon modes related to in- 
plane vibrations from Ca atom 9,10 . It has been found that R-3m phase with hexagonal ring transforms into layer 
orthorhombic Cmmm structure with five- and eight-membered rings at 18 GPa 11 . Further, the pressure dependence 
of T c studied in ac susceptibility measurements indicated that following an initial increase at the rate +0.39(1) K/ 
GPa, T c drops abruptly from 15Kto4KatlO GPa and no superconducting transition is observed above 2 K 
between 18 and 32 GPa 8 . The repressed superconductivity at about 32 GPa should be correlated to further structural 
modification, signaling the vicinity of new phase. So far, however, available structural information is scare, which 
partly hinders theoretical or experimental studies on highly compressed CaC 6 , since it is indispensable to discern 
stable structures so as to explore hidden physical property of CaC 6 at higher pressure. 

Results 

In this study, ab initio evolutionary algorithm (EA, USPEX code) 12 ~ 14 , which has been used to successfully predict 
high pressure phases of well-known materials, from element solids 15-16 to binary 17,18 or ternary 19 compounds, is 
chosen to address low enthalpy phases of compressed CaC 6 . Under compression, a novel high pressure phase with 
nanoporous characteristic, controlling a wide pressure range, presents superconducting feature with considerable 
T c of 14.7 K at 78 GPa. The re-emergence of superconductivity has been observed in various systems. Applying 
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magnetic field, Levy et al. 20 observed superconductivity in the ferro- 
magnet URhGe. Yuan et al. 21 observed two disconnected supercon- 
ducting phases in the pressure-temperature phase diagram of 
partially germanium-substituted CeCu 2 Si 2 . Very recently, Sun et 
al. 22 reported that an unexpected superconducting phase re-emerges 
in the superconducting iron chalcogenides above 1 1.5 GPa. Here, re- 
entrance of superconductivity in compressed CaC 6 provides new 
example to discern peculiar physical behavior of material at extreme 
conditions. 

Pressure induces phase transitions from hexagonal to orthorhom- 
bic structures. Three orthorhombic high pressure new phases were 
found. The enthalpies vs. pressure curves of stable structures are 
presented in Fig. 1 a). We checked stability of CaC 6 compared 
to its element solids using calcium's stable structures reported 



experimentally and theoretically 23 . No decomposition is observed 
in our calculations. The R-3m structure (referred as phase I) is the 
most stable structure below 19 GPa, followed by an orthorhombic 
C222 structure (referred as phase II, 19-39 GPa). The Pmmn phase 
(referred as phase III), a good metal (see supplementary Fig. SI), 
predominates in the pressure range of 39 to 126 GPa. Diamond-like 
carbon structure is adhered to the phase IV (Cmcm) for compressed 
CaC 6 above 126 GPa, showing a bad metallic feature (see supple- 
mentary Fig. S2) due to existence of most of saturated sp 3 hybridized 
carbons. Phase II with buckled five- or eight-membered ring's plane 
has very close enthalpy as of Cmmm structure with flat five- or eight- 
membered ring's plane reported by Csanyi et al. 11 . For these two 
structures, we cannot say which structure is more stable than other 
in term of total energy due to their very small structural difference. As 
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Figure 1 | (a) Relative enthalpy per chemical formula unit of compressed CaC 6 as a function of pressure, referenced to the Pmmn structure, (b) Structural 
phase transitional diagram under cold pressure, (c) Atomic arrangement of the R-3m, C222, Pmmn, and Cmcm structures of CaC 6 . The large ball 
represents the Ca atom. There are one, two, four, and three in equivalent carbon atoms in the R-3m, C222, Pmmn, and Cmcm structures, respectively. 
Under pressure, six-membered ring re-arranges into five- and eight-rings. Further, carbon polymerizes into three dimensional networks. Carbon mixing 
hybridization of sp 3 + sp 2 is observed in Pmmn and Cmcm phases. For Pmmn structure, sp 2 hybridized bonds are linked by sp 3 bonds, forming nanoporous 
structure, which is centered by calcium atoms chains. For Cmcm structure, the arrangement of carbon atoms is transformed into sp 3 framework structure 
(diamond-like behavior) with sp 2 bonds as borders and calcium ions have rearranged into 2D buckled plane. 
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pressure increases, buckling characteristics in phase II gradually dis- 
appears, making it transform into Cmmm symmetry. Phonon spec- 
trum calculations show that both C222 and Cmmm structures are 
dynamically stable. It is well known that non-hexagonal rings in a 
perfect graphene can be generated by introducing topological 
defects 24,25 . Here, external pressure and chemical pre-compression 
from calcium atoms provide the activation energy to drive di-vacan- 
cies defects formed, which leads to graphene layer in phase I trans- 
formed into octagonal and pentagonal carbon rings as observed in 
phase II. 

In phase III, there are four types of in-equivalent carbon atoms in 
the unit cell (referred as CI, C2, C3, and C4, see supplementary Table 
SI), constructing carbon nanoporous structure. Along y-axis dir- 
ection, carbon atoms form near quadrangular microporous, in which 
the size of porous lies in below 2 nm (about 0.61 nm and 0.59 nm at 
0 GPa and 39 GPa, respectively). CI and C2 (four carbon atoms in 
unit cell) show sp 3 hybridization, other two types of carbons C3 and 
C4 form zigzag chains in the middle of border with two different 
bonding lengths (1.426 and 1.467 A at 39 GPa, respectively, see 
supplementary Table S2), showing mostly sp 2 hybridization. Ca 
atoms stay at center of porous (coordination number 8, with C3 
and C4 atoms), forming ID atomic chain along y-axis direction. 
At phase transition pressure (39 GPa), the nearest neighbor distance 
between C and Ca atoms is about 2.31 A, which is slightly smaller 
than the value 2.489 A of Ca-Ca bonding length (i.e., lattice para- 
meter b), so there is stronger interaction between sp 2 hybridized 
carbon and calcium because of charge accumulation in the bond 
region (see supplementary Fig. S3). Strong covalent bonding between 
carbon atoms (sp 3 and sp 1 hybridization) bring about its high rigidity, 
thus having high hardness. To explore its rigidity, we calculate elastic 
constants (see Table 1) and phonon spectrum of Pmmn structure at 
zero pressure. No imaginary frequency (see supplementary Fig. S4) 
was observed together with elastic stability in conjuncture with 
elastic and dynamical stability of carbon host lattice (space group 
Cmcm) at zero pressure (see Table 1 and see supplementary Fig. S5), 
signaling its rigidity. We also observed a meta-stable high pressure 
phase (P2j/m) for compressed CaC 6 , in which carbon nanoporous 
structure includes armchair-type sp 2 carbon chains in the middle of 
border in porous. For phase IV, under external pressure and chem- 
ical pre-compression from Ca sub-lattice, carbon atoms integrate 
closely into diamond-type property's strip (a total of five layer carbon 
atoms), while Ca atoms is also rearranged as the bulked plane com- 
pared to centered sites of porous in phase III. In diamond-like strips, 
internal carbon atoms form distorted tetrahedrons with obvious sp 3 
hybridization (see supplementary Fig. S3), while surface carbon 
atoms mostly form sp 2 hybridization in armchair arrangement way 
along z-axis direction. The alternate carbon 3D framework and metal 
layers have an AtxB)3 arrangement, where the A and B represent the 
metal layers and a and fj the diamond-like strip layers. The formed 
covalent three dimensional quasi-layer network (Ca layer-diamond- 
like layer-Ca layer- diamond-like layer) in phase IV largely enhances 
its incompressibility and hardness in comparison with phase III. 
Using suggested formula by Simunek and Vackaf 26 , the estimated 
hardness is 34 GPa for Pmmn CaC 6 at zero pressure, 45 GPa for 
Pmmn CaC 6 at 39 GPa, and 91 GPa for Cmcm CaC 6 at 126 GPa. 



The atomic arrangements of competing structures are shown in 
Fig. 1 c). The nearest neighbor distance between carbon atoms is 
1.444 A for the R-3m phase at 0 GPa, 1.417 A for the C222 phase 
at 19 GPa, 1.426 A for the Pmmn phase at 39 GPa, and 1.460 A for 
Cmcm phase at 126 GPa, which are close to the atomic distances in 
graphite. Interestingly, we found that only small change is observed 
for these distances as pressure increase, indicating very stiff chemical 
bonding between carbon atoms. For phase III from zero pressure to 
39 GPa, carbon-carbon length only decreases 3 ~ 4%. The stiffness 
of C-C bonds corresponds to high vibration frequencies in phonon 
dispersions (discussed later). 

Carbon nanoporous structure is very important because of its 
potential applications in many fields, such as H 2 absorption, nano- 
devices, supercapacitors and lithium ion batteries, as well as chemical 
and bio-sensors. Structurally, the carbon foam, a three-dimensional 
porous including both sp 2 and sp 3 hybridized atoms, can be thought 
of as assembled from graphene planes interconnected rigidly with 
one another at a specific angle, forming a linear chain of sp 3 bonded 
atoms along the junction 27 . Previously, supposed hybrid carbon foam 
structures with a mixed sp 2 /sp 3 bonding character 28,29 have received 
much less attention for lack of direct experimental observation. 
Recently, Tomanek suggested carbon foam may form under non- 
equilibrium conditions near gain boundaries of a carbon-saturated 
metal substrate 30 . Here, (2, 2) zigzag carbon foam observed in phase 
III of CaC 6 points out new easy-realized way experimentally via 
powerful high pressure tool. Electronic structure calculation shows 
that (2, 2) zigzag carbon foam is a semiconductor with narrow band 
gap (—0.49 eV, see supplementary Fig. S6). In view of nanoporous 
scientific and technological importance, we in detail explore host 
lattice's structural behavior (carbon lattice in phase Pmmn) at zero 
pressure. Removing calcium from lattice, relaxed host carbon sub- 
lattice shows higher symmetry (parent group Cmcm), which attri- 
butes to CI and C2 as well as C3 and C4's symmetric degeneracy. For 
Ca guest sublattice, it forms body-centered tetragonal lattice with 
Immm symmetry. When adding calcium atoms into host lattice, 
interaction between Ca and C leads to lattice's symmetry decrease. 
Most exciting is that carbon nanoporous structure relaxed remains 
stable mechanically and dynamically at zero pressure (see Table 1 
and supplementary Fig. S5). At zero pressure, carbon foam equilib- 
rium structural parameters are listed in supplementary Table SI. The 
densities, carbon-carbon bond lengths, and carbon-carbon bond 
angles of carbon foam and phase III of CaC 6 as well as experimental 
values of graphite and diamond are listed in supplementary Table S2. 
The foam contains two-thirds 3-fold coordinated carbon atoms, 
labeled sp 2 , and one-third 4-fold coordinated carbon atoms, labeled 
sp 3 . The mass density of the optimized foam structure, p = 2.589 
g/cm 3 , lies in-between the experimental values for graphite and dia- 
mond (see supplementary Table S2), which is highest density among 
reported carbon foams so far. Our carbon foam has peculiar struc- 
tural behavior. There are three types of covalent bonds, sp 2 -sp 2 , sp 2 - 
sp 3 , and sp 3 -sp 3 and five types of bond angles, sp 2 -sp 2 -sp 2 , sp 2 -sp 3 -sp 2 , 
sp 2 -sp 3 -sp 3 , sp 3 -sp 2 -sp 3 , and sp 3 -sp 3 -sp 3 , as presented in Table supple- 
mentary S2, which leads to non-hexagonal cross-section at junction. 
From supplementary Table S2, it can be seen that sp 2 -sp 2 bond length 
in carbon foam is shorter than that in graphite, again indicating its 



Table 1 | Elastic constants, bulk, shear, and Young's modulli (unit: all in GPa) of stable phases of CaQ as well as carbon foam (labeled Q) at 
zero pressure 
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stronger rigidity, which differs from previous reported carbon 
foam 27,29 . The size of nanoporous is often measured by its border 
length. For (2, 2) zigzag carbon foam observed here, it has 3.925 A 
of width of tube-wall. The pore diameters (diagonal lengths) of nano- 
porous are 5.935 A and 5.408 A. Further, nanoporous's size can be 
modified by change numbers of carbon six-membered rings of bor- 
ders, which helps to tune specific surface area and band gap so as to 
facilitate requisite applications in purification and separation, 
adsorption, and catalysis 31 . This shape modification of nanoporous 
has been observed in other compressed GICs KC 8 and LiC 6 . Using Li 
as guest, we obtained (2, 2) nanoporous and triangular nanoporous, 
and using K, two kinds of nanoporous discovered (see supplement- 
ary Fig. S7). We can expect that more nanoporous solids can be 
exposed by other metal guest (such as Sc, Ba, Cs, and Yb et al). 
Structure decides property, so further studies remains desirable to 
uncover structural evolution behavior of GICs under compression. 
Also we found that some carbon nanoporous systems with peculiar 
shape observed here present negative linear compressibility 31 . To add 
calcium atoms into (2, 2) zigzag carbon foam brings about porous 
distortion due to strong hybridization between calcium d-electrons 
and sp 2 hybridized carbon p-electrons (see supplementary Fig. SI). 
As is shown in supplementary Table S2, there are five bond lengths 
and nine bond angles for carbon atoms. Carbon-carbon bond is 
slightly increased in length because of calcium atomic insertion but 
presents strong rigidity. 

The small density of states (DOS) value at Fermi level for phase II 
exactly explains disappearance of T c of CaC 6 at the pressure range 
from 18 to 32 GPa experimentally 8 . Phase III has high DOS value at 
Fermi level compared with that of phase I (see supplementary Fig. S8) 
stimulates our intense interest on its potential superconductivity, 
since the electronic DOS at the Fermi level sensitively controls T c . 
Energy bands of phase III at 39 GPa were shown in Fig. 2 (a). There 
are four energy bands, labeled as 1, 2, 3, and 4, crossing the Fermi 
level along several high symmetry directions in BZ. The most inter- 
esting feature is the appearance of a small region (TY direction in BZ) 
of one two-fold degenerating flat band lying just at the Fermi level 
around the Y point. The occurrence of flat and steep slopes near the 
Fermi level, which is a favorable condition for enhancing Cooper pair 
formation, has been thought to be essential to superconductivity 32,33 . 
The calculated DOS shows that a strong hybridization between the 



Ca-<i electrons and sp 2 carbon-p electrons (C3 and C4-p electrons) 
appears in the valence band from —2.78 to 0 eV and in the conduc- 
tion band from 0 to 4 eV. Electrons from sp 3 carbon (CI and C2)-p 
electrons construct sp 3 chemical bonding, lying mostly valence band 
from— 16 eVto— 2.78 eV. Dominant contributions to the DOS near 
the Fermi level come from the Ca-d electrons and sp 2 carbon-p elec- 
trons (C3 and C4), whereas the Ca-p, Cl-p, and C2-p electrons make 
minor contributions to the DOS at the Fermi level (see supplement- 
ary Figure SI). Figure 2 (b) shows the Fermi surface (FS) of the Pmmn 
structure at 39 GPa. The labeled numbers (1-4) correspond to the 
labeled energy bands shown in Fig. 2 (a). The FS consists of two 
electron-hole oblate tube with an '8'Tike section along the S-Y dir- 
ection (all four bands), two hole tubes (bands 3 and 4) near the S-X- 
U-R plane, and two hole pockets (band 2, T-Z direction). Obvious 
near parallel pieces of the Fermi surface exists, which favors to 
enhance the electron-phonon coupling. Both the Fermi pockets 
and the obvious Fermi nesting create strong electron-phonon inter- 
actions in CaC 6 . By further analyzing the wave functions of the four 
bands at the high symmetry points, we conclude that electron-like 
Fermi surfaces arise mainly from the Ca-c/ z 2 and Ca-d x 2 states, 
whereas hole-like Fermi surfaces attributes to the p x and p z states 
from C3 and C4. 

Discussion 

Phonon dispersions and the partial phonon density of states 
(PPHDOS) of phase III are shown in Fig. 3. The flat bands observed 
along S-X-U-R direction may serve as evidence of its local two-dimen- 
sional feature. Many soften modes are observed along the T-Z-T-Y 
and T-R directions (see also supplementary Fig. S9). As seen in the 
PPHDOS (in the right panel of Fig. 3), the calcium atomic vibrations 
mostly contribute to the phonon dispersions below 385 crrT 1 . The sp 3 
carbon (CI and C2) atomic vibrations mainly dominate narrow fre- 
quency range around 1200 cm -1 . The phonon dispersions between 
385 and 889 cm 1 as well as above 1226 cm 1 are mainly due to the 
contributions of the sp 2 carbon (C3 and C4) atomic vibrations. 

Using the phonon linewidth 7q , we can identify the contribution 
to the electron-phonon interaction parameter X from each mode 

(X 0 ) based on the relation — ; r -> ■ Here, MO) is the density 

of states at Fermi surface. We notice that large values of Xq lie at the 



1600 




Figure 3 | Phonon dispersion along the high-symmetry directions of the Brillouin zone (left panel) and the partial phonon density of states (PPHDOS) 
(right panel) of the Pmmn CaC 6 at 39 GPa. 
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Figure 4 | (a) The Eliashberg phonon spectral function a 2 F{oj) and 
electron-phonon integral X(co) of the Pmmn CaC 6 structure at 39 GPa. 
(b) Superconducting transition temperature versus pressure for CaC 6 . 

high symmetry points T, S, X, U, and R. For instance, we found that 
at the R point, ^(co) originates mainly from contributions of the 
vibration modes at the frequency ranges of 242 to 305 cm" 1 . The 
calculated spectral function a 2 F(co) and integrated X(io) of Pmmn 
CaC 6 at 39 GPa are plotted in Fig. 4 a). The vibrations below 
385 cm 1 provide the major contribution to X (about 0.45). The 
low-frequency phonons (below 385 cm -1 ), which mostly involve 
the Ca and sp 2 carbon (C3 and C4) atoms yield 54.5% of the total 
X value. The phonons between 385 and 889 cm 1 mainly from sp 2 
carbon atoms, contribute to 32.4% of total X. These results indicate 
that the calcium and sp 2 carbon atoms in the Pmmn structure dom- 
inate superconductivity in CaC 6 , due to the prominent contributions 
to the electron-phonon interaction. Phonons from the Ca and sp 2 
carbon atoms together with the electrons from the Ca-d and sp 2 C- 
p states provide the strong electron-phonon coupling necessary for 
strong superconductivity in Pmmn CaC6. 

The Allen and Dynes modified formula 34 was used to estimate the 
superconducting transition temperature T c from the value of X deter- 
mined above. Taking a typical value of 0.115 for /(* along with the 
calculated coi og of 642 cm -1 , we obtained a T c of 3.48 K for Pmmn 
CaC 6 at 39 GPa. Applying pressure to Pmmn CaC 6 clearly leads to X 
first enhanced remarkably, then weaken and thus T c corresponding 
responded (see supplementary Table S3). The calculated highest 
value of T c is up to considerate 14.7 K (at 78 GPa). At 117 GPa, T c 
is only 4.74 K. The calculated decrease of T c agrees with occurrence 
of another new phase (Cmcm). The change of T c with increasing 
pressure in CaC 6 is summarized in Fig. 4 b). Phase III has consistent 
superconducting mechanism with phase I. Re-emerged supercon- 
ductivity in CaC 6 is closely related to structural transformation. 



Carbon nanofoam centered calcium atom chains in phase III brings 
about good metallic feature and strong electron-phonon coupling 
and thus superconductivity. Pressure-induced atomic hierarchical 
arrangement in phase IV, presenting diamond-like carbon strips 
(mostly sp 3 carbon) and bucked calcium atomic layers, decreases 
its conductivity and thus represses superconductivity. By tuning 
electronic order (orbital or spin), pressure or magnetic field-driven 
the re-emergence of superconductivity has been observed in various 
systems including ferromagnet URhGe 20 , heavy fermion compound 
CeCu 2 Si 2 21 , and iron chalcogenides 22 . 

To conclude, for compressed CaC 6 , carbon atom's arrangement is 
transformed from preceding two-dimensional (2D) layer into three- 
dimensional (3D) framework, while calcium ions themselves from 
3D network turn into 2D buckled plane. A novel high pressure phase 
(Pmmn, two molecules/cell) with nanofoam carbon shows super- 
conducting behavior due to enhanced electronic distribution at 
Fermi level as well as strong electron-phonon coupling, capturing 
critical superconducting temperature of 14.7 K at 78 GPa. The pres- 
sure-induced re-emergence of superconductivity in CaC 6 signals a 
rich physics for graphite intercalation compounds under high com- 
pression. Although compressed LiC 6 holds nanoporous carbon in its 
stable high pressure phases but it does not show the superconducting 
behavior. Therefore, further theoretical and experimental efforts are 
well worth to discern physical nature of GICs under compression. 

Methods 

Stable structures of GICs under cold compression were searched using evolutionary 
algorithm in combined with VASP code 35 based on density functional theory within 
the generalised gradient approximation with the exchange-correlation functional of 
Perdew Burke Ernzerhof 36 , employing the projected augmented wave (PAW) 3738 
method where 2s 2 2p 2 and 3s 2 3p 6 4s 2 are treated as valence electrons for C and Ca 
atoms, respectively. For the crystal structure searches, we used a plane- wave basis set 
cutoff of 700 eV and performed the Brillouin zone integrations using a coarse k-point 
grid. The most interesting structures were further relaxed at a higher level of accuracy 
with a basis set cutoff of 1000 eV and a k-point grid of spacing 2n X 0.018 A -1 . 
Iterative relaxation of atomic positions was stopped when all forces were smaller than 
0.001 eV/A. The lattice dynamic and superconducting properties of Pmmn CaC 6 
were calculated by the Quantum ESPRESSO package 39 using Vanderbilt-type ultra- 
soft potentials with cutoff energies of 55 Ry and 500 Ry for the wave functions and the 
charge density, respectively. The electronic Brillouin zone (BZ) integration in the 
phonon calculation was based on a 12 X 16 X 12ofMonkhorst-Packk-pointmeshes. 
The dynamic matrix was computed based on a 2 X 4 X 2 mesh of phonon wave 
vectors. The electron-phonon coupling was convergent with a finer grid of 48 X 64 X 
48 k points and a Gaussian smearing of 0.01 Ry. 
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